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Abstract We measured mobility patterns that describe walking trajectories 
of individual Me'Phaa peasants searching and collecting fuelwood in the forests 
of "La Montana de Guerrero" in Mexico. These one-day excursions typically 
follow a mixed pattern of nearly-constant steps when individuals displace from 
their homes towards potential collecting sites and a mixed pattern of steps of 
different lengths when actually searching for fallen wood in the forest. Dis- 
placements in the searching phase seem not to be compatible with Levy nights 
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described by power-laws with optimal scaling exponents. These findings how- 
ever can be interpreted in the light of deterministic searching on heavily de- 
graded landscapes where the interaction of the individuals with their scarce 
environment produces alternative searching strategies than the expected Levy 
flights. These results have important implications for future management and 
restoration of degraded forests and the improvement of the ecological services 
they may provide to their inhabitants. 

Keywords Levy flights • deterministic walks • Me'Phaa • fuelwood • ecological 
restoration • human mobility 



Introduction 

Mobility patterns of humans searching for resources se em to follow specific 



statis tics known to physicists as anomalous diffusion (Klaflcr an d Sokolov . 



2005) . From the prehist oric tribes in Europe searching for quality stone guar 



rels |Brantinghaml . 120061) to the hunter-gatherers of Africa ( Brown et all . 12007 



2010T 1. it is becoming clear that distances traversed when humans travel, mi- 



grate, disperse and explore a territory is not quite a Gaussian s t ochastic pro- 
cess but one described by scale- free statistic s (JBrockmann et all . l2006t iFlores , 
l2007t iGonzalez et~al 120081: iRhee et at l201l[ ). 

Understanding the non-Gaussian mobility nature of humans is very im- 
portant for realistically explaining diverse patterns in humans such as mating 
(gene flow), cultural changes, foraging, spread of disease and migration pro- 
cesses, all of which have potential and unexpected impacts on the daily life of 
individuals and the complexity of societies. Mobility drives the dynamics of 
encounters among humans and influences the way vital resources and ecosys- 
tem services are used, for example food and energy sources. In this study, we 
measured the patterns of displacements of individual peasants when searching 
and collecting fuelwood in the forests of "La Montana de Guerrero" in the 
southern Mexican state of Guerrero. We discuss evidence that these patterns 
do not conform to Levy flight statistics dominated by power-law probability 
distributions of the form P(l) ~ l^^ with scaling exponent fi (1 < fi < 3) 
and I being the length of the displacements. The absence of Levy flight pat- 
terns is not, however, evidence against Levy foraging as a human mobility 
phenomenon but rather that, in this case, it is the result of foragers inter- 
acting with an environment characterized by the scarcity typical of degraded 
landscapes. 

Conventional optimal foraging theory states that animals and humans 
search for resources guided by decision s made in order to maximize a given 
quan tity, usually the net energy gain ( Charnovl . Il976t IStephens and Krebd . 



1986). The efficiency of the search is then in proportion to the behavioral 



repertoire of the foragers who may employ, to varying degrees, cognitive tools 
such as memory and information processing when facing the searching space. 
Most searching models fall in two main groups: (1) random searching where 
the forager does not have previous information of the prey field such as the 
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spatial distribution or density and (2) deterministic searching where the for- 
ager is able to use memory maps of the territory and therefore has relatively 
good previous information about the prey field richness and location. 



Random searching 



Random searching is by no means e xclusively perfo rmed by simple low-profile 
organisms (Matth aus et all . 1201 lb lHavs et all . 12011 ). it may be used even by 



humans under specific environment constraints, for example when searching 
for a plane that has crashed in the open ocean at an unknown location. It is an 
optimal strategy when there is no previous information about where the target 
or prey is located. Recently, optimal random se arching theory has been devel- 
oped further around the so called Levy flight s (|Shlesinger and Klafterl . l2000t 
Chechkin et al . 2008; Viswanathan et all . 120111) . It has been demonstrated that 
step-lengths generated by a power-law probability distribution is an optimal 
searching strategy when the scaling exponent is fi ~ 1 when the prey is re- 
moved (destructive scenario) o r \x = 2 when t he prey is not removed (non- 
destructive scenario) (jViswanathan et all . Il999i) . 

In the destructive scenario there is no reason to re-visit previous sites 
such that this searching strategy is conveniently done by "ballistic motion" 
where successive step-lengths are very large in size and are only separated 
by a random change of angle. In practice, it corresponds to the strategy of 
"follow a straight line until you find something" . Ballistic searching was first 
suggested to exist in the very long rectilinear flights of the wanderi ng albatross 
but w as later called into question (jEdwards et all . 120071 but see also lBover et al 
(2008)). A non-destructive searching scenario, in which the prey can regenerate 
or reappear in the same location at a later time, is optimally performed when 
the forager alternates a few long steps with more abundant short steps. This 
scenario with /i ~ 2 was found to exist in many biological examples but has 
also been rec e ntly r eviewed and questioned citing methodological inadequacies 
([James et all 120111 ). Since rece nt unquest i oned examples o f Levy flights are 
continuously being published ([Sims et all . l2008t Irlavs et all . 120111 ). the state 
of the art of this research area has attracted much interest and is rapidly 
changing. 



Deterministic searching 



Animals and particularly humans may rely on cognitive tools and memory 
maps in order to perform more efficient searches on a territory previously 
known or similar to one previously known. A common strategy is "go to the 
nearest richest place" in which the forager makes a decision on where to go 
by weighing the profit offered by the site against the energy invested in reach- 
ing it, depending on the forager's present location. A well-studied problem 
related to this deterministic searching mode is the tourist problem in which a 
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tourist wishes to visit a number of cities without repeating them and where 
some cities ar e more attractive than othe r s, posing a challenging optimiza- 
tion problem ([Stanley and BuldvrevL l200ll llima et all . l200ll : ICampiteli et al 
20061 ). Scale- free movement statistics are known to be an emergent property in 
models of decision-making intelli gent agents performing deterministic searches 
in heterogeneous enviro nments (jBover and Larraldd . l2005t I B over et all . l2009t 
Bover and Walshl . l2010i ). Emergent scale- free mobility patterns have been sug- 
gested to exist in animals ranging from mic rozooplankton to spider monkeys, 

among others (|Revnoldd . l2008t iBover et all . |2004 IgOO ^I). 

In the case of the spider monkeys search patterns ([Ramos-Fernandez et all . 



20041 ) it has been suggested that these arise as a consequence of the interaction 
of the searching individuals with the environment in such a way that charac- 
teristics l ike density and sp atial richness distribution influence the searching 
strategy ( Bover et all . 120061) . Different searching options may be displayed in 



different quality habitats. In rich habitats, resources such as food are so abun- 
dant at close ranges that long-distance displacements are not expected. At 
the other extreme, when the environment is poor with scarce resources, there 
is no point in traveling far away because the energy invested in long trips is 
not justified since the gain will be the same that when moving locally. Again 
long-distance displacements are not expected. The situation is very different 
when the environment is neither exceptionally rich nor poor . In this case Levy 
like d isplacement statistics are robustly expected to emerge ( Bover et all2006l 
20091) . 

Levy-like statistics have unexpected consequences for social foragers. For 
example, in spider monkeys, it has been shown that complex social struc- 
tures spontaneously emerge from the interactions of individuals moving under 
such conditions. These s ocial networking even have a small-world dynamics 
( Ramos- Fernandez et all . I2006T) . In the case of the Dobe Ju/hoansi hunter- 
gatherers mobility patterns, it was suggested that the Levy foraging influ- 
enced the way tribes visited and remained for a time in specific foraging sites. 
When individuals from different tribes came into contact in these sites, cul- 
tural exchanges and inter-tribe couple mating (gene flow) took place. In this 
way a complex soc ial structure emerged paired to the local use of resources 
(JBrownet all . l2007l) . 



The Me'Phaa of La Montana de Guerrero 



Me'Phaa is a pre-Columbian indigenous ethnic group that inhabits the re- 
gion known as "La Montana de Guerrero" (MG) in the state of Guerrero 
in Southwestern Mexico. Me'Phaa towns scatter the area around the city of 
Ayutla de los Libres (16°54' N, 99° 13' W). The Me'Phaa were until very re- 
cently knew as "Tlapanecas" which is pejorative Nahuatl term from the time 
the Aztec Empire ruled the region. MG is one of the poorest and least de- 
velo ped areas in Mexico, with Human Develop ment Index (HDI) of around 
3.2 ( Morales- Hernandezl 12006c iTaniguchiL 120111 ). comparable to some areas of 
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sub-Saharan Africa despite its close proximity to the world-class resort of Aca- 
pulco. A typical Me'Phaa individual has no access to health services, schools, 
paved roads, telecommunications, or electricity, and their situation has been 
worsened by the recent inciden ce of militarization, social conflict, and violence 
(JHebertl . |200J ICamachol 120071 ). 



Fuelwood usage and availability 

The landscape of MG is hilly with forest patches mostly on hilltops (see Fig- 
ure [1]). Seasonally tropical dry forest is the most prevalent type of vegetation, 
but oak and pine forests are present in the highlands (around 1890 msl). The 
environment is sev erely degraded with high erosion and deforestation rates 
(JLanda et all . Il997l) . MG is known to have been subject to o c cupat ion and 



anthropogenic impact for as long as 2000 years (jBerrfo et all . 120061 ) . which 



has caused the present pattern of degradation. Fuelwood is mostly composed 
of species belon ging to the gene r a Aca cia, Leucaena, Lysiloma, Prosopis and 
Pithecellobium ( Cervantes et all . 119981) . Quercus and Pinus species are also 



collected when the searching is done in the highlands. 

Many Me'Phaa peasants do not have any access to modern energy sources 
and most continue to rely primarily on fuel wood for low-tech residential uses 
such as cooking and heating. MG is regarded as one of the main hot spots of 
fuelwood usage in Mexico and so requires priori ty attention since thi s practice 
continues to stress the already degraded forest (jGhilardi et all 120071 ). 



Material and methods 

Twelve voluntary Me'Phaa peasants were trained to operate GPS data loggers 
in order to register positions along a trajectory in the field when searching 
for fallen fuelwood. The learning process posed no difficulty and after few 
minutes of manipulating the devices the volunteers showed adequate skill when 
operating them. Specific instructions were given to switch on the devices when 
leaving home on a search and collection excursion and to switch them off when 
returning home. The devices were configured to log data automatically as soon 
as the first accurate position fix was obtained. The GPS data loggers were 
worn on a neck lanyard and protected by a plastic bag from dust and rain. 
All volunteers were males in their 20s to 50s and the searching for fuelwood 
was done alone with no supervision or accompaniment by the research team. 
The peasants were left on their own to carry out the searches whenever and 
wherever they do so normally. Only fallen wood was collected, and no trees 
were felled. Normally, fallen wood is picked up and carried by hand until a 
sufficient amount is collected and tied together into a bundle carried over the 
shoulders. All volunteers belonged to NGO for organic agriculture production 
based in the city of Ayutla de los Libres (Xuajin MePhaa AC). Participants 
live in scattered communities in the mountains surrounding the city, separated 
from each other by several kilometers. 
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We used Holux-M241 data loggers (Holux Technology Inc.) equipped with 
a MKT GPS-chipset capable of storing 100,000 points. The devices were set to 
record successive points at intervals of 10 seconds. All devices were recovered 
and the data transferred to a computer for analysis. The data was analyzed 
with an ad-hoc program we developed that implemented a Haversine algo- 
rithm to calculate distances between two geo-referenced points on earth. The 
participants recorded a total of 114 field trips, but only 10 of them contained 
search-related data. These recordings contained a total aggregated amount of 
3386 displacements, from which 1231 displacements corresponded to searching 
displacements with a maximum value of 28.4 meters. All displacements less 
than 1 meter were ignored when analyzing the displacements since such dis- 
tance in 10 seconds is approximate to a waiting interval. These values however 
were considered when analyzing the waiting times. 

All statist ical analysis was conducted using the R Open Source Statisti- 
cal Language (|R Development Core Teaml . 120091 ) on a Linux machine (Ubuntu 
11.04). Step-lengths and waiting-times dis t ributions were fitted w ith two alter- 
native statistical models (JNewmanl . l2005t Edwards et all . 120071 ): a power-law 
l~v- and the exponential e~ xl . Maximum Likelihood Estimation (MLE) was 
used for evaluating model parameters and loglikehood values. In a set x of 
random variables of size n: x = {xi, X2, ■ ■ ■ , x n }, the scaling exponent /i of the 
power-law is giving by: 



jj, = 1 + n 



L 4 t 



(1) 



In the case of the exponential distribution, the A parameter value is given 



by: 



A 



E 

.1=1 



(2) 



In both cases x a is the minimal value of x from which the probabil- 
ity distribution hold; in our case x a was fixed as the minimun x value, so 
that the entire time-series was evaluated. A Model Selection approach involv- 
ing an Akaike Information Criterion (AIC) was conducted in order to iden- 
tify among the two stati s tical models the one that better explain the data 
([Burnham and Anderson . I2002I ). The AIC of model i is defined as AICi — 
—2Li + 2Di, where L is the log-likchood of the fit and D is the number of 
parameters of the model. The model with minimal value of AIC is regarded as 
the most parsimonious and is normally considered the most likely to explain 
the data ( Burnham and Andersonl 120021) . 



Results 



The one-day search and collecting excursions typically follow a mixed pattern 
of nearly-constant steps when individuals displace from their homes towards 
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the forest (a "ballistic phase" ) and a pattern of steps with a fat-tail like distri- 
bution when actually doing the searching for fallen wood (see Figured]). This 
pattern of mobility behavior is very common in nature and has received severa l 
name s (JKnoppien and Reddingiusl . Il985t lLomholt et all l2008t iBenichou et al , 
20111 ). It is characterized by an adaptive switching between searching behav- 
iors depending on the availability and location of the targets. The first phase is 
dominated by a simple displacement of the forager to the area of interest and 
a second phase follows when a potentially resource-rich location is reached. In 
this second phase, active and careful search behavior is initiated. We ignored 
the displacements pertaining to the first "ballistic" transportation phase and 
analyzed only the displacements in the second active search phase. 

Step-lengths distribution was fitted with two alternative statistical mod- 
els using Maximum Likelihood Estimation, a power-law and the exponential. 
The parameter values found are giving in Table [1] (see Figure [3]). A Model 
Selection approach points to the exponential model as the most parsimonious, 
having the lower AIC value. On this basis, we conclude that the Me'Phaa 
searching displacements are essentially a Brownian-like motion. On the other 
hand, the waiting times, that is the time intervals with no displacements, were 
analyzed with the same procedure as above (see Figure |4] and Table [I]). In 
this case the lower AIC value suggests that the most likely model explain- 
ing the data is the power-law model (fi = 2.6), meaning that the resource is 
scattered in patches with amounts of fuelwood following a scale-free probabil- 
ity distribution, a result that is in agreement with the known fact that the 
size distribution of forest t rees and their associated biomass is scale-invariant 
(|Enauist and Niklasl . I2O0I . 



Discussion and summary 



Mobility is increasingly becoming a problem of interest for anthropology ( Brown et all . 



20071 : iRicherson and Bovdl . 12008). Of particular interest is how individual for- 



aging movements influence the pattern of social contact s that in turn deter mine 
the emergence of comp lex social structures, in human s (jBrown et all . I2007T ) and 
non- human primates ( Ramos-Fernandez et all . [2006J). The patterns of social 
contacts are also important for understanding th e spatial dynamics of socia l 
phenomen a such a s cultural and social changes dRicherson and Bovdl . 120081) . 
gene flow (jSlatkinl . Il973l; ISokal et at Il989t IWakelevl. Il999l). disease spread- 
ing a nd vaccination strategies ( Miramontes and Luauel 120021 : lMao and Bian. 
120101 ). among others. Human mobility and the environment are also entan- 
gled in a complex matrix of feedback interactions where mobility directly im- 
pacts the availability of resources and specific ecosystem services at various 
scales and where environmental fluctuations and landscape degradation may 
also negatively influence the pa tterns of human mobility and social behavior , 
even e nhancin g social conflicts (JHomer- Dixonl . Il994r iRaleigh and Urdall . 12007 ; 
iBurke et all . I2009J: iHsiang et all l201ll ). 
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For many years, it was thought that human migrations and individual mo- 
bility could be described by Gaussian probability distributions; however recent 
studies suggest that human mobility may be better explained by anomalous 
diffusion where the statistics of displacements follow power-law distributions in 
the form of Levy flights. The origins of Levy statistics in human displacement 
is therefore an issue of increasing interest. Due to the fact that humans search 
using cognitive tools for decision-making, most of their foraging behavior is 
aimed at optimizing a cost /benefit ratio, as conventional foraging theory pre- 



dicts (jCharnovl . 119761 ) . Such optimization is at the core of a modern approach 
that argues for the existence of a deterministic behavior of searching where the 
spatial distributions of the target field (richness and density ) cause the sponta- 



neous emergence of the Levy fo raging in intelligent agents (pantos et all . 12007 : 



Bover , 2008; Bover et all . I2009T) . When the environment is scarce in resources, 



deterministic searching force non-Levy patterns of displacements because the 
Levy stable distributions start converg ing into a Gaus s ian one when the scal- 



ing e xponent fj, of the power law is > 3 (jChechkin et all . 12008c iNurzaman et al 



20111 ). Such non-Levy patterns are the best local strategy for such environ- 
ments but these are non-global optimal solutions. 

It is estimated that nearly 2.5 billion people in developing countries world- 



wide make use of fuelwood in order to meet their residential energy needs ( IE At 



2009). Most of these people live in areas subject to strong environmental pres 
sure. This is the case of the Mexican Me'Phaa peasants in "La Montana de 
Guerrero" in Mexico. This impoverished indigenous group has inhabited the 
area since pre-Columbian times and so have made use of the ecosystem ser- 
vices since then with non-persistent large-scale management strategies. We 
have studied the mobility patterns of the Me'Phaa peasants when searching 
and collecting fuelwood in the field. We have found a mixture of foraging be- 
haviors that consist of a long trip (up to several kilometers on foot) from their 
homes to the collecting sites. This phase is composed of nearly constant dis- 
tance steps that mostly follow the path of roads or paths. This is similar to 
a "ballistic behavior" where there is no active searching at all. When arriv- 
ing to an interesting site that contains a significant amount of fuelwood, the 
movement behavior is replaced by an active search composed of a mixture of 
abundant small steps al ternated with few long steps in a f ashion reminiscent of 
intermittent searchin g (JKnoppien and Reddingiua . Il985c lLomholt et all l2008t 



Benichou et all 120111 ). The statistical distribution of the step lengths is ex- 
plained better by an exponential model. 

What is the meaning and the origin of an apparently non-optimal Brownian- 
like movement pattern in the Me'Phaa searching process? First we should ex- 
amine in more detail the searching behavior as performed by the individuals. 
When arriving in the area to be searched, peasants do have a fairly good view 
of where the fallen wood lies around since the density of trees is low in a de- 
graded seasonally tropical dry forest (also, there are no other major physical 
obstacles). This means that the individuals would not move around search- 
ing randomly as if blind. Instead, they will move to where they see (tens of 
meters) there is fallen wood of good size and in good amount. This searching 
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Table 1 MLE 


parameter values of two models 


Model 


Parameter 


Loglikc 


AIC 


Step-lengths 


Power-law 
Exponential 


fi= 1.97 
A = 0.03 


-3336.00 
-3274.37 


6674.00 
6550.75 


Waiting-times 


Power- law 
Exponential 


,u = 2.66 
A = 0.47 


-162.44 
-261.01 


326.89 

524.03 



behavior is then repeated and it matches the pattern of a deterministic search 
as explained in the introduction. We rule out the behavior leading to a Levy 
distribution with op timal scaling exponents p, — 1 or fi = 2, typical of random 
searching scenarios ([Viswanathan et all 120111) . Instead we emphasize that, in 
models of deterministic searching, it has been argued already that environ- 
ments with scarcity of resources may lead to Brownian-like displacements. This 
does not mean that the search is inefficient but that a Brownian-li ke searching 
pattern bec omes locally optimal under such extreme conditions ( Bover et all 
20061 120091) . 

Bro wnian-like behavior is also known to o ccur in other non- human pri- 
mates (jSchreier and Grovel l2010l ISueurl 120111 ). stressing the importance of 
considering the role of the environment in influencing the mobility patterns 
of foragers. Spatial distribution, density an d resource abundance can not be 
separ ated from the study of forager mobility (|de Jager et all 1201 lllBover et al , 
2011) , espec i ally w hen these involve decision-making and optimization efforts 
(JBoveret alll2009l) . 

Fuelwood search and collection by the Mexican Me'Phaa peasants is char- 
acterized by large amounts of energy invested in traversing large distances 
walking in order to profit poorly, since the collection is a low-tech activity 
limited by the carrying capacity of the individuals. Therefore the searching 
behavior of the Me'Phaa despite of involving elaborate decision-making and 
complex interactions with the environment is far from a global optimal solu- 
tion. This situation may be reversed or improved if Me'Phaa peasants could 
actively modify the scarce nature of their forests by means of adopting forest 
restoration practices that would increase the availability of wood close to their 
homes. Research on the restoration of ecosystem services in the region is an 
activity in progress following this study. 
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Fig. 1 Panoramic view of the landscape of "La Montana de Guerrero" in Mexico, one of 
the poorest regions in the country. The image at the left shows the hilly nature of the 
seasonally tropical dry forest currently characterized by large rates of deforestation and soil 
erosion. The image shows one example of a walking trajectory (T) performed by a GPS- 
equipped peasant when traveling to collect fuclwood. The area of large white spots at the 
bottom of the image is a village (V) of Mc'Phaa peasants of about 20 houses (17.1388969 
N, -98.9851532 W, 1314 msl). At the end of the recorded trajectory there is a shift in the 
searching behavior that becomes an area-restricted active search. The behavior is shown in 
the square window that has been enlarged in the image at the right. This last image shows 
in detail, the actual displacements over the terrain when searching. Images are an overlap 
of the GPS positions on Google Earth™ imagery. Images courtesy of Google Inc. 









(d) 




Fig. 2 Two examples of travel and search trajectories as recorded every ten seconds by 
walking Mc'Phaa peasants (a and c). These images show in detail the two phases of the 
searching and collecting of fuelwood. The first phase is a simple "ballistic" displacement 
towards the potential collecting site and may span up to several kilometers. When the 
peasants arrive to the site of interest, the behavior is shifted as shown in the two enlarged 
figures (b and d). The pattern of mobility is then replaced by an entangled succession of 
many short steps and few large steps. 
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Fig. 3 Statistics of the step lengths. This is a log-log plot of the histogram of the distances 
traveled by the peasants versus the number of such distances, recorded every ten seconds. 
The straight line shows the power-law fit as estimated using the MLE method. The power-law 
has an scaling exponent fi = 1.97. The curved line is the exponential tested as an alternative 
model over the same interval (see Table [lj- The exponential model better explains the data 
as suggested by an AIC model selection. 
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Fig. 4 Statistics of waiting times. This is the log-log plot of the histogram of the waiting 
times versus their frequency, recorded every ten seconds. The straight line shows the power- 
law fit as estimated using the MLE method. The power-law has an scaling exponent fi = 2.66. 
The curved line is the exponential tested as an alternative model over the same interval (see 
Table [Q. The AIC model selection suggests that the power-law is the most adequate model 
to explain the data. 
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